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ABSTRACT 

The effect of halo substructures on galaxy rotation curves is investigated in this paper using a simple 
model of dark matter clustering. A dark matter halo density profile is developed based only on the 
scale free nature of clustering that leads to a statistically self-similar distribution of the substructures 
at galactic scale. Semi-analytical method is used to derive rotation curves for such a clumpy dark 
matter density profile. It is found that the halo substructures significantly affect the galaxy velocity 
field. Based on the fractal geometry of the halo, this self-consistent model predicts an NFW-like 
rotation curve and a scale free power spectrum of the rotation velocity fluctuations. 
Subject headings: dark matter — galaxies: dwarf — galaxies: general — galaxies: ISM — galaxies: 
structure — ISM: kinematics and dynamics 
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1. INTRODUCTION 

Observatio nal and theoretical s tudies of galaxy rota- 
tion curves (|Rubin fc Ford! 1 19701) . velo c ity d ispersions 
of elliptical galaxies dFaber fc Jacksonl 1 10T6f) , barvon 
fracti on in clusters (jClowe et al.l 20061 : IMassev et al.l 
20071) . gravitational lensing, structure formation, CMB 
power spectrum etc. indicate that a significant frac- 
tion of the mas s of the Universe h as no electromag- 
netic interaction (|Bertone et al.| [2005). At galactic scale, 
the main evidence for the existence of such gravitat- 
ing matter with very high mass to light ratio comes 
from observations of galaxy rotation curves. The 
gravitational potential derived from the observed ro- 
tation curve can not be explained only by the vis- 
ible mass with any reas o nable mass to light ratio 
(IRubin fc Ford [19701 [19851 ISohi^[l99l ISofue fc Rubir] 
20011: ISpano et al.ll2008D . Though there are mod ified the- 
ories of the gravitation and other suggestions dMilgrom l 
19831: iBekenstein fc Milgroml 119841: iSanderi Il986t iFahrl 
1990; Sanders 19971; iBrownstein fc Moffat! 120061) to ex- 



plain this anomaly, the dark matter concept is widely 
accepted to be a simpler explanation of these observa- 
tions. Though observations of galaxy rotation curves es- 
tablished the existence of the dark matter at galactic 
scale, there is no general agreement on the nature and 
various properties (like mas s distribution) of this major 
const ituent of the Universe (|Navarro et al.|[l996t Ide Blokl 
2005). Potentially, galaxy rotation curves can also shed 
light on some aspects of the dark matter properties and 
thus help us in deriving a better understanding of the 
nature of the dark matter. One such key aspect is the 
density profile or the mass distribution of the dark mat- 
ter at galactic scale. 

On the theoretical front, there are various mass dis- 
tribution models of the dark matter halo. The isother- 
mal profile, the Navarro-Frenk- White (NFW) profile and 
some variant of these two profiles are used extensively in 
both theoretical and observational studies. An almost 
flat observed rotation curve outside the core of a galaxy 
led to the dark m atter halo model with 1/r 2 isother - 
mal density profile ( Beecma n et al.lfT991b iBurk crt 1995). 
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The density at the centre is not finite in this model. 
However, it is possible to impose appropriate boundary 
conditions to derive a non-singular solution with den- 
sity p(t)nis = Po r c/( r c + r2 )' where po is the central 
density and r c is the "core radius". At large radius, 
this profile will be similar to 1/r 2 singular isothermal 
profile. A cut-off radius r max is also required to be im- 
posed to keep the total mass of the halo finite. How- 
ever, the non-singular isothermal profile was not very 
successful in explaining the observed rotation curves. A 
significant progress in the field was the introduction of 
an alternative density profile, the NFW density profile, 
p(t)nfw — Po r l/[ r ( r c + r ) 2 }, where p c and r c are the 
charac teristic density and characteristic ra dius respec- 
tively H^EoeEaLl EUl, [1992 UEBWmI)- 

The NFW profile has been quite successful in explain- 
ing the observed rotation curves for many galaxies. Over 
the time, a class of variant of the isothermal and the 
NFW profiles are also introduced t o fine tune the agree- 
ment of theory and observation s (|Burkertl [l995; Zhao 
119961 : iFukushige fc Makinoll2001[ ). All these density pro- 
files are of more or less similar merit, as far as explaining 
the rotation curve is concerned. However, detailed anal- 
ysis shows that the NFW profile has a number of prob- 
lems. For example, observations of gal axies and clus- 
ters suggest a centra l flat density core (Ide Blokl 120071 : 
iKuzio de Narav et al.1 120081). whereas, in NFW profile, 
dark matter density has a central cusp with a logarith- 
mic slope sa —1. Numerical simulations show that the 
presence o f baryon can change the mass distrib ution sig- 
nificantly ([El-Zant et al.ll2004t [Shlosmanll2010l and ref- 
erences therein), and for isothermal cusp, minor merger 
and dynamic al friction may lead to a shallower central 
density slope (|Romano-Dfaz et al.ll2008l) . The NFW pro- 
file also seems not to fit the observed rotation curves of 
the dark matter dominated low sur face brightnes s galax - 
ies and low mass dwarf galaxies (|de Blokl 12001 I2005T) . 
A more serious issue with the NFW profile is that the 
profile is derived by fitting analytical function to the 
density distribution derived from numerical simulation 
of structure formation with dark matter. This, in a 
sense, lacks a proper physical understanding. Though 
some physical insight of the NFW profile have recently 
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come from further high resolution dark matter simula- 
tions (e.g.. lTavlor fe Navarroll200l . the issue is far from 
being settled. Moreover, even if these cosmological sim- 
ulations are very successful on large scale, at galactic 
scale there are u nsolved issues, like the "angular momen- 
tum problem" (iSommer-Larsen. Gelato fe Vedell 119991: 
Sommer-Larsen fe Dolgovl 120011: iBurkert fe D'Onghia 



2004) or the "missing s atellite problem" (|Klypin et al.l 



1999; Mo ore et al.lll999fl . which are yet to be addressed 



Here, we have developed a self-consistent model of the 
dark matter halo substructure distribution at galactic 
scale to explain observed NFW-like rotation curves. Sup- 
port for the existence of dark matter subs tructures has 
mainly come from numerical simulations dGiocoli et al 



20081: iMadau et al.ll2008l:lSpringel et al.ll2008l : lElahi et al 
2009a HLudlow et al.ll2009h . But, there are strong obser 
vational hints lik e flux anomalies and time delays in grav- 
itational lensing (lChenll2009t iKeeton fe Moustakasll20"ol 



iVegetti fc Koopmansll2009l; IXu et al.ll2009fl. or enh anced 
gamma rays and le ptonic cosmic rays (|Elahi eFaIi r2009b: 
iPinzke et al.l 120091 ). indicating the presence of substruc- 
tures. The present model is based on the assumption 
of a scale free nature of the dark matter clustering that 
leads to a statistically self-similar distribution of the halo 
substructures at galactic scale. It is shown that a sim- 
ple fractal model of the dark matter halo substructure 
predicts an NFW-like rotation curve. Such a model also 
predicts a scale free power spectrum of the rotation ve- 
locity fluctuations. The model is described in Section [21 
and the results are presented in Section |3l Possible limi- 
tations of our analysis are discussed in Section0] Finally, 
we summarize and present our conclusions in Section [5] 

2. THE MODEL 

2.1. Assumptions, Parameters and Constraints 

Assuming that the dark matter clustering has a scale 
free nature (i.e., there exist halo substructures of a wide 
range of mass), the density profile can b e descri bed as 
a combination of a smooth radial profile p(r/r c ) and a 
stochastic part Sp. Here r c is a characteristic "core" 
radius and ~p contains information of the density vari- 
ation at scales larger than or comparable to r c . For the 
purpose of this work, we have used a simpler model of 
this density distribution which, however, retains all rel- 
evant key features. In this simplified model, we have 
assumed that each dark matter clump has a number of 
smaller clumps around it. Each of these smaller clumps 
are in turn just a scaled down version with even smaller 
clumps around them. As a result, the whole structure 
has an approximate spherical symmetry and a statistical 
self-similarity. It is assumed here that all these clumps 
have non-singular isothermal density profile where the 
central density (po), the core radius (Y c ), the cut-off ra- 
dius (r max ) and the halo to subhalo distance (D) are 
scaled down accordingly. However, this specific density 
profile is not a crucial assumption in our model and the 
individual substructures may have any non-singular den- 
sity profile p(r/r c ), where r c is some characteristic ra- 
dius. Essentially, this is a fractal structure with three 
parameters: (i) n, the number of small clumps around 
any clump, (ii) f r , spatial scaling factor for core radius, 
cut-off radius and distance and (iii) f p , central density 
scaling factor between any clump and its next smallest 



clumps. A fractal is a fragmented and irregular geometri- 
cal shape wi th exact or stocha stic self-similar structures 
at all scales ((Ma ndelbrot 1983). Independent of the value 
of f r and n, the (Hausdorff) fractal dimension of such a 
structure is 3 for nff < 1. However, since at each iter- 
ation, the linear size of clumps scales by a factor f r and 
the mass scales by a factor n, the local mass dimension 
for any substructure level is D m = —log{n)/log(f r ) over 
a certain range of scales. A mass dimension of D m for 
a medium implies that the mass enclosed in a sphere of 
radius r in such a medium will be M(r) — kr Dm . So, 
for the N th substructure level, M^{r) — k^r ™ over a 
range of scales depending on N, n, f r , D and r max . Note 
that this range is different for different substructure lev- 
els. Hence, for the complete structure, the total mass 
M(r), which is the sum of Mjv(r) of all the substructure 
levels, will not have a simple power law radial depen- 
dence. But, the dark matter halo mass function will still 
be a power law, N(m) oc m~ Q , where the power law 
index a = —log(n)/log(ff,f p ) is a more physically mo- 
tivated parameter of this model and can be constrained 
from theoretical and numerical analysis of dark matte r 
structure formation (jGao et al.ll2004t iZemp et al.l l2009). 
The density distribution can be written as 
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p(r) = Pb e + 2^ 2^ Ps(po,i,rc,i,r m ax,,' 

i=0 j=l 



n,j) (1) 



where pbg 

Ps(.P0,i: Tc.ii Tmax.i - ' /. / ; 

ual substructure centred at r; 



is background density and 
,ri_j) is density profile of individ- 
Ei t j with central density 
{po.i), core radius (r Cj i) and cut-off radius (r maXy i). 
Considering the self-similarity of this model, 



Po,i — fp x Pcu-i 
^c,i — fr x r c i 
^max.i — fr X r ma 



n,j = ?Vi,fc + (f r x A-i x a) 



(2) 



where k — n , d is a unit length vector with 
random orientation and the initial set of parameters 
Po,o,r C:(h r maX:0 = p ,r c ,r max is for the largest subhalo 
centered at the origin. In principle, pb e may be a smooth 
function of r. But, since we are assuming it to be a small 
background density threshold, its effect on the final ro- 
tation curve is not very significant. So, for simplicity, we 
have assumed pbg to be constant over the radius of our 
interest. This structure is shown schematically (without 
any randomness for the shake of clarity) in the left panel 
of Figure [TJ After introducing randomness in angular 
position of the subhalos, one realization of such a struc- 
ture with n = 7 and /,. = 0.33 is shown with two and 
four substructure levels in the middle and right panel re- 
spectively. This model can be considered as a simplified 
representation of the scale free, clumpy density structure 
of dark matter above a small threshold density at the 
galactic scale. 

We note that the parameters for this model are con- 
strained to a good extent by various physical consider- 
ations. Assuming that the structure is extended down 
to the infinitely small scale, to avoid a divergence of 
the total mass, the quantity nfffp must be less than 
unity. Similarly, for a halo with cut-off radius r max , halo 
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Figure 1. Left: Schematic representation of the halo substructure with three subhalos around any big halo. Note that both the cut-off 
radius (r max ) and the halo to subhalo distance (D) are scaled down by a factor f r for each substructure level. Middle and Flight Halo 
substructure with seven subhalos around any big halo and with f r = 0.33. The middle and right panels show the structure with two and 
four substructure levels respectively. 



to subhalo distance D must be greater than or equal 
to (1 — fr)r ma x/0- ~ 2/ r ) to avoid any overlap. Also 
note that the density of any subhalo at the cutoff radius 
p{Tmax) scales as f p . So, we adopt f p = f.O for further 
analysis to ensure a constant density at the cutoff radius 
for all the substructures. In this simplified model, we 
have assumed that all the parameters like D, po, r max 
etc. are identical for all the subhalos for a particular sub- 
structure level and n, f r and f p remain constant for all 
the levels. In a realistic situation, however, all these pa- 
rameters may have some random variation. As a result 
of these fluctuations, the halo mass function, the radial 
mass distribution, and in turn, the rotation curves, are 
expected to be somewhat smoother than that are pre- 
dicted from this analysis. 

2.2. Tidal stability 

A much stronger constraint on the parameters comes 
from the consideration of the stability of this structures 
preventing tidal disruption by invoking self-gravity. Con- 
sidering a rigid object of mass m and radius r at a dis- 
tance d from a bigger object of mass M and radius R, 
from the standard Roche limit consideration, D should 
be ~ r(2M/m) 1 / 3 to avoid tidal disruption of the smaller 
body. In the case of any halo and its immediate subhalo, 
m/M = f$ and r/R = f r imply d « 1.25i?. For even 
smaller subhalo structures with r/R — /*, the mass is 
scaled accordingly m/M = / 3fe and d m 1.25R assures 
stability. Hence, a distance 

D > 1.25a; 1 ~ ^ r max (3) 
l zj r 

will make the whole structure stable. Here x is a fudge 
factor and we use the value of a; = 1.1 to accommodate 
non-rigid density clumps. For a given value of D, the 
number of substructure n is also constrained to be 
4^D 2 _4(l-/ r ) 2 



n < 



n(f r D/l - f r f 



to avoid any possible overlap of subhalo structures. 
2.3. Virial stability 



A detailed virial stability analysis requires numerical 
simulation of the dynamics of such a density distribution 
to get the time-averaged dynamical quantities. But, a 
simple ensemble average virial scaling analysis may be 
used to constrain the central density po for a set of model 
parameter. As the whole structure is assumed to have 
an approximate spherical symmetry, average kinetic and 
potential energies, (T) and (V), for thin spherical shells 
of radius r and thickness Sr will be 



(T) = -map M 



-47rr 2 <5rp(r)oD M 



= _ GM(r)m = _ ^ (5) 

r 

where ctdm is dark matter velocity dispersion, M(r) is 
total mass within radius r and v c (r) = {GM /r) 1 ' 2 is the 
scale dependent virial velocity equivalent of the "circular 
velocity" for rotating disk. Since the rotation curve has 
a roughly constant value fo at large radius, the ratio 
2(T)/\(V) | will tend to the equilibrium value of 1 at large 
radius if ctdm ~ ^o- Using the minimum value for D from 
equation ([3]), the maximum extent of the structure Rmax 
will be £)/(! — f r ) and the average density will be 



(P) 



3(l-2/ r ) 3 (s -tan^s) 
(1.25xs) 3 (l - nff) 



-Pc 



(6) 



where s = r max /r c . Note, however, that this is a frac- 
tal structure with significant porosity. So, the average 
density of any individual clump is higher than (p) by a 
factor of (1 -nf3)(R max /r max ) 3 . Now, for global stabil- 
ity of the whole structure, R m ax should be less than or 
equal to the radius within which the virial equilibrium 
is maintained. Average density within this virial radius 
r v i r or r2oo should be » 200 times more than the critical 
density p cr u — 3H 2 /8ttG. For a choice of model param- 
eters n, f r , s and x, this will constrain the lower limit of 
the central density po so that (p) > 200 p c nt- For indi- 
vidual substructures, both mass and volume scale as / 3 , 
keeping the average density constant. This implies that 
stability for one substructure level ensures stability for 
all other levels. 
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Figure 2. Predicted rotation curve for the fractal substructure 
model (black) and the best fit NFW profile to that (red). This 
is for n = 35, f r = 0.25, f p = f.O and r m ax = 6.0r c . Rotation 
curves for non-singular isothermal sphere with and without a cut- 
off (green and blue curves respectively) as well as for NFW halo 
(magenta curve) are also shown for a comparison. Radial distance 
and rotation velocity are scaled by r c and vq (rotation velocity at 
the furthest radial distance) respectively. 

3. RESULTS 
3.1. Rotation curve 

Since the halo density distribution is significantly 
clumpy, the velocity field for such a system is also ex- 
pected to have fluctuations at all scales. However, due to 
approximate spherical symmetry of the clump distribu- 
tion, the average rotation velocity over a spherical shell 
at radius r will still be (v c (r)) ~ (GMr/r) 1 / 2 , where M r 
is the total mass within this radius. As pointed out in 
subsection 12.31 using the virial stability argument, the 
virial velocity or, equivalently, the "rotation" velocity is 
expected to be approximately same as the local velocity 
dispersion. This derived rotation curve for the fractal 
model is found to be NFW-like at large radial distance. 
At small radius, by construction, the rotation curve is ob- 
viously exactly same as that of a non-singular isothermal 
halo. This is shown in Figure[5J Radial distance is scaled 
by r c and rotation velocity is scaled by vq, which is the 
rotation velocity at the maximum radius plotted. The 
black curve is the predicted rotation curve for the fractal 
substructure model with n — 35, f r — 0.25, f p = 1.0 
and r max = 6.0r c and the red line is the best fit NFW 
profile to that. The background density threshold pb g is 
assumed to be negligible in this model. Rotation curves 
for non-singular isothermal sphere with and without a 
cutoff (green and blue curves respectively) as well as for 
NFW halo (magenta curve) are also shown in Figure [2] 
for a comparison. The velocity scaling ensures that the 
total mass encompassed by the maximum radius is same 
for all models. 

The effect of the background density threshold pb g is 
shown in Figure [3] Here, the rotation curve is evaluated 
for different pb g keeping all other parameters same as in 
Figure [5J For scaling, we have used the density p(r max ) 
at the cutoff radius r max = 6.0r c . Different curves in 
Figure Hare for p bg /p(r max ) = 0.00, 0.03, 0.10 and 0.30 
(black, red, green and blue curve respectively). The gen- 
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Figure 3. Predicted rotation curve for the fractal substructure 
model for different background density threshold pb g . All other 
parameters are same as in Figure [2] Different curves are for 
Pbgl p{ r max) = 0.00, 0.03, 0.10 and 0.30 (black, red, green and 
blue curve respectively). 



eral NFW-like nature of the rotation curve and the radial 
fluctuations remain unchanged. But depending on the 
value of pb g , rotation curve at large radius may be rising, 
flat or declining. Note that the derived rotation curves 
shown in Figure [3] are with the simple model of a con- 
stant pb g ■ In reality, pb g is expected to be decreasing with 
increasing r, giving rise to a rotation curve somewhat in- 
termediate between the extremes shown in Figure [3] 

3.2. Velocity fluctuations 

The predicted rotation curve due to the clustering of 
dark matter subhalos is very similar to the observed rota- 
tion curve and the empirical NFW rotation curve. How- 
ever, unlike the NFW model with a smooth radial den- 
sity distribution, the present model predicts a significant 
fluctuation of rotation velocity in both angular and ra- 
dial directions. Since the underlying density field, which 
gives rise to this velocity fluctuations, is scale free, the 
velocity fluctuation power spectrum is expected to be a 
power law. Though the fractal model has many free pa- 
rameters, the only relevant parameters for the index of 
this power law are n and f r while the rest of them will 
just introduce different multiplicative scaling. This pre- 
diction can be easily verified from high spatial and spec- 
tral resolution observation of neutral hydrogen of normal 
galaxies. Note that part of this fluctuations will cancel 
out for the spherically averaged rotation curve and hence 
it is important to use the full velocity field to search for 
such scale free fluctuations. It is also important to note 
that fluctuations of the velocity field of the hydrogen gas 
will have contributions from the local density perturba- 
tions of the disk. But the scale dependence of these per- 
turbations may makes it possible to disentangle the scale 
free perturbations due to the halo substructures. There 
is some indication of such a power law scaling of the 
rotation velocity fluctuation power spectra from direct 
observation and analysi s of the H I 21 cm velocity field 
of some nearby galaxy (|Dutta et al.ll2010l ). We leave a 
more detailed treatment of this aspect to a future work. 
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4. DISCUSSION 

Throughout this analysis we have assumed that the 
clump distribution in this model has an approximate 
spherical symmetry. This assumption is most likely to 
be untrue in reality. Due to small count, bigger sub- 
structures are more likely to cause departure from spher- 
ical symmetry. In this situation, rotation velocity field 
will also have strong azimuthal asymmetry. Observa- 
tionally, about half of the spiral galaxies show some de- 
gree of kinematical l opsidedness (|Richter fc Sancisilll994l: 
Havnes et aT1ll998t ISwaters et al.l l2008t iSofue fc Rubin! 
20011 : Uog fc Combesfl 2009) which may originate from 
tidal distortion partly due to deviation from spherical 
symmetry of the underlying g ravitational potential of the 
lopsided dark m atter halo (jChakrabarti fc Blitd 120091 : 
ISaha et aD I2008H . Due to their large number, smaller 
substructures are likely to have relatively more symmet- 
rical distribution. So, the velocity fluctuations spectrum 
is expected to show a power law scaling at large spa- 
tial frequency (large fc, i.e., small physical scale), and a 
departure from power law at small k. 

We have also assumed in this analysis that the baryons 
will not significantly affect the galactic dynamics. How- 
ever, some recent numerical studies (e .g. [Weinber g et al.l 
120081: iRomano-Di'az et al.l [20081 I2009D have shown that 
the presence of baryons have effects like flattening the 
central cusp, reducing the halo triaxiality, and intro- 
ducing bias in clustering with environment density. Ef- 
fectively, baryon dissipation may destroy the similar- 
ity between different sca les of clustering. However, 
I Romano- Diaz et al.l (|2010f) have found that phenomenon 
like "efficient feedback from stellar evolution and the cen- 
tral supermassive black holes" may counte rbalance this 
effect to some extent. On the other hand, iKnebe et al.1 
(|20Tot) have found no significant effect of baryonic physics 
on properties like shape and radial alignment of substruc- 
tures. Numerical simulations also suggest that, even in 
the presence of baryons, the subhalo mass function re- 
mains a power law tho ugh the power law index changes 
from -0.99 to -1.13 (jRomano-Diaz et al.l 12010^ . It is, 
hence, important to keep in mind that, in a realistic sit- 
uation, baryon dissipation may alter the velocity fluctu- 
ations causing a significant departure from the scale free 
velocity fluctuations. 

We note that the key result of this analysis, that is 
an NFW-like rotation curve for the fractal model, is 
not crucially dependent on the exact density profile of 
individual halos. A variety of density profile without 
any central singularity (variant of non-singular isother- 
mal sphere) will leads to a similar NFW-like rotation 
curve. This strongly suggests that the clustering proper- 
ties of the dark matter particles dominantly govern the 
radial density profile of the halo. Finally, these assembly 
of substructures will give rise to a flat NFW-like rota- 
tion curve for normal galaxies. But, for more dark mat- 
ter dominated low mass dwarf galaxies and low surface 
brightness galaxies, depending on the exact form of the 
non-singular density profile, the dominant contribution 
of the central big halo may make the rotation curve for 
such galaxies intrinsically different which is c onsistent 
with observational results (|de Blokll2002[ I2005D . 

5. CONCLUSIONS 



1. The dark matter halo substructures at galactic 
scale is found to significantly affect the rotation 
curve. 



2. A self-consistent model of statistically self-similar 
hierarchical dark matter substructures predicts an 
NFW-like rotation curve, though each clump has 
a non-singular isothermal density profile. This 
NFW-like rotation curve emerges out of the frac- 
tal geometry and is independent of specific density 
profile of individual clumps. 

3. The model predicts a scale free power spectrum 
of the rotation velocity fluctuations which can be 
observationally verified. 

4. The model also provides some plausible explana- 
tion of the observed intrinsic difference between the 
dark matter halo density profile of normal galaxies 
and dark matter dominated low mass dwarf galax- 
ies and low surface brightness galaxies. 



We are grateful to Rajaram Nityananda for many use- 
ful comments on an earlier version of this paper. We 
thank Susmita Chakravorty, Aritra Basu, Sanjay Bhat- 
nagar, Abhirup Datta, Prasun Dutta, Sanhita Joshi and 
Chandreyee Sengupta for helpful discussions. We are 
also grateful to the anonymous referee for useful com- 
ments and for prompting us into substantially improv- 
ing this paper. NR is a Jansky Fellow of the National 
Radio Astronomy Observatory (NRAO). This research 
was supported by the NRAO. The NRAO is a facility of 
the National Science Foundation (NSF) operated under 
cooperative agreement by Associated Universities, Inc. 
(AUI). 

Disclaimer. This is an author-created, un-copyedited 
version of an article accepted for publication in the As- 
trophysical Journal. IOP Publishing Ltd is not respon- 
sible for any errors or omissions in this version of the 
manuscript or any version derived from it. The definitive 
publisher authenticated version will be available online at 
http : / / iopscience . iop . org/. 



REFERENCES 



Begeman K. G., Broeils A. H., Sanders R. H., 1991, MNRAS, 
249 523 

Beken'stein J., Milgrom M., 1984, ApJ, 286, 7 
Bertone G., Hooper D., Silk J., 2005, Physics Reports, 405, 279 
Brownstein J. R., Moffat J. W., 2006, MNRAS, 367, 527 
Burkert A., 1995, ApJ, 447, L25 

Burkert A. M., D'Onghia E., 2004, Penetrating Bars through 
Masks of Cosmic Dust: the Hubble Tuning Fork Strikes a New 
Note, ed. D. L. Block et al. (Dordrecht: Kluwer), 341 

Chakrabarti S., Blitz L., 2009, MNRAS, 399, L118 

Chen J., 2009, A&A, 498, 49 

Clowe D., Bradafi M., Gonzalez A. H. et al., 2006, ApJ, 648, L109 
do Blok W. J. G., 2002, A&A, 385, 816 
de Blok W. J. G., 2005, ApJ, 634, 227 

de Blok W. J. G., 2007, in Island Universe, ed. R. S. de Jong 

(Dordrecht: Springer), 89 
Dutta P., Roy N., Majumdar S., 2010 (in preparation) 
Elahi P. J., Thacker R. J., Widrow L. M., Scannapieco E., 2009a, 

MNRAS, 395, 1950 
Elahi P. J., Widrow L. M., Thacker R. J., 2009b. Physical Review 

D, 80, 123513 

El-Zant A., Hoffman Y., Primack J., Combes F., Shlosman I., 

2004, ApJ, 607, L75 
Faber S. M., Jackson R. E., 1976, ApJ, 204, 668 
Fahr H. J., 1990, A&A, 236, 86 



6 



Fukushigc T., Makino J., 2001, ApJ, 557, 533 

Gao L., White S. D. M., Jenkins A., Stoehr F., Springcl V., 2004, 

MNRAS, 355, 819 
Giocoli C, Tormcn G., van den Bosch F. C, 2008, MNRAS, 386, 

2135 

Hayncs M. P., Hogg D. E., Maddalcna R. J., Roberts M. S., van 

Zee L., 1998, AJ, 115, 62 
Jing Y. P., 2000, ApJ, 535, 30 
Jog C J., Combes F., 2009, Phys. Rep., 471, 75 
Keeton C R., Moustakas L. A., 2009, ApJ, 699, 1720 
Klypin A., Kravtsov A. V., Valcnzuela O., Prada F., ApJ, 1999, 

522, 82 

Knebe A., Libeskind N. I., Knollmann S. R., Yepes G., Gottlobcr 

S., Hoffman Y., 2010, MNRAS, 405, 1119 
Kuzio de Naray R., McGaugh S. S., de Blok W. J. G, 2008, ApJ, 

676, 920 

Ludlow A. D., Navarro J. F., Springel V., Jenkins A., Frcnk C S., 

Helmi A., 2009, ApJ, 692, 931 
Madau P., Dicmand J., Kuhlen M., 2008, ApJ, 679, 1260 
Mandelbrot B., 1983, The Fractal Geometry of Nature, W.H. 

Freeman, New York 
Massey R., Rhodes J., Ellis R. ct al., 2007, Nature, 445, 286 
Milgrom M., 1983, ApJ, 270, 365 

Moore B., Ghigna, S., Governato F. et al., 1999, ApJ, 524, L19 
Navarro J. F., Frcnk C. S., White S. D. M., 1996, ApJ, 462, 563 
Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 
Pinzke A., Pfrommer C, Bergstrom L., 2009, Physical Review 

Letters, 103, 181302 
Richtcr O.-G., Sancisi R., 1994, A&A, 290, L19 



Romano-Diaz E., Shlosman I., Heller C. H., Hoffman Y., 2008, 
ApJ, 685, L105 

Romano-Diaz E., Shlosman I., Heller C. H., Hoffman Y., 2009, 
ApJ, 702, 1250 

Romano-Diaz E., Shlosman I., Hoffman Y., Heller C. H., 2010, 

ApJ, 716, 1095 
Rubin V. C, Ford Jr. W. K., 1970, ApJ, 159, 379 
Rubin V. C, Ford Jr. W. K., 1980, ApJ, 238, 471 
Saha K., Levine E. S., Jog C. J., Blitz L., 2009, ApJ, 697, 2015 
Sanders R. H., 1986, MNRAS, 223, 539 
Sanders R. H., 1997, ApJ, 480, 492 

Shlosman I., 2010, Galaxy Evolution: Emerging Insights and 

Future Challenges (ASP Conf. Ser. 419), ed. S. Jogee et al., 39 
Sofuc Y., 1996, ApJ, 458, 120 
Sofue Y., Rubin V., 2001, ARA&A, 39, 137 
Sommer-Larsen J., Dolgov A., 2001, ApJ, 551, 608 
Sommcr-Larscn J., Gclato S., Vedel H., 1999, ApJ, 519, 501 
Spano M., Marcelin M., Amram P. et al., 2008, MNRAS, 383, 297 
Springcl V., Wang J., Vogelsbcrgcr M., Ludlow A., Jenkins A., 

Helmi A., Navarro J. F., Frenk C. S., White S. D. M., 2008, 

MNRAS, 391, 1685 
Swatcrs R. A., Schocnmakcrs R. H. M., Sancisi R., van Albada T. 

S., 1999, MNRAS, 304, 330 
Taylor J. E., Navarro J. F., 2001, ApJ, 563, 483 
Vcgctti S., Koopmans L. V. E., 2009, MNRAS, 400, 1583 
Weinberg D. H., Colombi S., Dave R., Katz N., 2008, ApJ, 678, 6 
Xu D. D., Mao S., Wang J., Springel V., Gao L., White S. D. M., 

Frcnk C. S., Jenkins A., Li G., Navarro J. F., 2009, MNRAS, 

398 1235 

Zemp'M., Diemand J., Kuhlen M. et al., 2009, MNRAS, 394, 641 
Zhao H., 1996, MNRAS, 278, 488 



